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ABSTRACT 

The aeroelastlc stability of a circulation control rotor blade undergoing three 
degrees of motion (flap, lag, and torsion) is investigated in forward flight. Quasi- 
steady strip theory is used to evaluate the aerodynamic forces; and the airfoil charac- 
teristics are from data tables. The propulsive and the auxiliary power trims are 
calculated from vehicle and rotor equilibrium equations through the numerical integra- 
tion of element forces in azimuth as well as in radial directions. The nonlinear time 
dependent periodic blade response is calculated using an Iterative procedure based on 
Floquet theory. The periodic perturbation equations are solved for stability using 
Floquet transition matrix theory. The effects of several parameters on blade stability 
are examined, including advance ratio, collective pitch, thrust level, shaft tilt, 
structural stiffnesses variation, and propulsive and auxiliary power trims. 


INTRODUCTION 

The airfoil on a circulation control (CC) rotor typically has quasl-elliptlc profile 
and uses a tangential wall jet ejected over the rounded trailing edge to produce cir- 
culatory lift; (fig. 1). Due to the Coanda effect, the air remains attached at the 
rounded trailing edge and the stagnation point shifts to the lower surface. The lift of 
a CC airfoil can be controlled by jet momentum as well as by geometric Incidence. It 
is possible to achieve high lift coefficients (four to five) with CC airfoils. Also, 
the aerodynamic center due to blowing circulation is near half-chord. For a general 
bibliography of circulation control see Englar and Applegate (1984). The application of 
CC technology to full-scale rotor design is currently being investigated (Linden and 
Biggers, 1985). The cyclic lift control of a CC rotor is obtained by a cyclic modula- 
£> tlon of blowing. This eliminates the need for cyclic pitch and results in a simplified 
hub design. Collective lift control is obtained either by collective blowing or by 
collective pitch; see figure 2. Important features of a CC rotor are high thrust capa- 
bility at reduced tip speeds and easy implementation of a higher harmonic control 
system. One area of concern, however, is the effect of CC aerodynamics on blade dyna- 
mics. 

Johnson (1985) documented recent developments in the dynamics of advanced rotor 
systems. Few attempts have sbeen made to examine the aeroelastlc stability of a CCR 
rotor blade (Chopra and Johnson, 1979; Chopra, 1984-85). An aeroelastlc stability ana- 
lysis of a CC rotor blade in hover was conducted by Chopra and Johnson in 1979. Three 
degrees of motion were considered: rigid flap, lag, and feather rotations about hinges 
at the blade root. The CC airfoil characteristics were represented In terms of analyti- 
cal expressions. It was shown that the trailing edge blowing can have an Important 
Influence on blade dynamics and must, therefore, be addressed in rotor design. Chopra 
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(1984) also examined the aeroelastic stability of flap bending, lead-lag bending, and 
torsion of a CC rotor hingeless blade in hover using a finite element formulation. The 
CC airfoil characteristics were taken from data tables. Again, the stability results of 
hingeless rotors showed that the blowing has an important influence on blade dynamics. 
The finite element formulation was extended to analyze the aeroelastic stability of a 
bearingless rotor blade in hover (Chopra, 1985). 

The objective of the present work is to examine the aeroelastic stability of a CC 
rotor blade in forward flight. For this, a simple flap- lag-tors ion blade model con- 
sisting of three degrees of motion is considered. Quasi-steady strip theory is used to 
obtain aerodynamic forces. The effect of unsteady aerodynamics is introduced approxima- 
tely through dynamic inflow modeling. The effects of pneumodynamics (Watkins et al. , 
1985) and centrifugal pumping in the pressure duct are included to calculate the jet 
momentum at a radial station. 

The propulsive trim is calculated iteratively from the vehicle nonlinear equilibrium 
equations. Three force equations (vertical, longitudinal, and lateral) and two moment 
equations (pitch and roll) are obtained by numerically integrating the element forces 
both along the azimuth as well as the radial directions. The trim solution gives the 
rotor control setting and the vehicle orientation for a prescribed flight condition. The 
blade steady response is then calculated from nonlinear periodic blade equations using 
an iterative procedure based on Floquet theory. For stability, the blade motion is 
assumed to be a small perturbation about the steady response, and the linearized 
periodic blade equations are solved using Floquet transition matrix theory (Panda and 
Chopra, 1985). Stability results are calculated for typical CC rotor blades for several 
flight conditions. 


FORMULATION 

The blade is assumed to undergo three degrees of motion: rigid flap, lag, and 
feather rotations about hinges at the blade root, with hinge springs to obtain the 
desired natural frequencies. The hinge sequence is flap inboard, followed by lag, and 
then feather outboard. The flap angle 6 is positive up, the lag angle £ is positive aft 
(opposite to rotation), and the feather angle 0 is positive leading-edge up. The struc- 
tural equations of motion include the inertial forces about three hinges and are docu- 
mented by Panda and Chopra (1985). In general, terms up to second-order are retained in 
the flap and lag equations and terms up to third order are retained in the feather 
equation. 

The aerodynamic forces are obtained using quasisteady strip theory. The section 
lift, drag, and moment about the mid-chord (per unit span) are 

L - j pV 2 cC Jl (a,C i| ) 

D - I pV 2 cC d (a.C^) (1) 
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The aerodynamic coefficients C^, C^, and C m# 5 are from data tables; the numerical 
values for these coefficients are available at small steps, Aa of 3 degrees and AC^ of 
1/200. These coefficients depend on the airfoil geometry, including slot height, and 
are also a function of angle of attack a, blowing momentum coefficient C , and local 
Mach number. However, in the present work, the effect of compressibility (Mach number) 
is neglected. The is defined as 


C 


li 
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is the jet momentum, 
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For an incompressible flow, using an isentropic 
momentum coefficient can be related to the local 
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duct pressure P^. 
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where h/c is slot height-to-chord ratio (typically 0.002) and is duct gage 

pressure. The dynamic pressure at a radial station is 

q - - p (OR) 2 (- + vsin<|0 2 
2 R 

For a compressible flow with a subsonic jet condition, the blowing momentum coef- 
ficient (Rogers et al«, 1985) is 


s " 2 £ ( v m j 2 for (M j < L) 


(4) 


where is the jet Mach number and M^ is free-stream Mach number. These relations, 
eqs. (3) and (4), for incompressible and compressible flows are valid only for the 
unchoked flow condition (when P./p^ is less than 1.892. For a P^/P^ larger than 1.892, 
the flow becomes choked in the nozzle (slot). Then, the jet momentum coefficient is 
obtained as 
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For both cases of unchoked and choked flows, the jet Mach number can be calculated in 
terms of duct pressure 


/ p *r 

M - /s [(li) T - 1 ] (6) 

J r oo 

For an accurate representation of the blade internal duct pressure characteristics 
at an arbitrary local rotor disc location, pneumodynamic considerations are included in 
the analysis (Watkins et al., 1985). Between the duct pressure at a radial station and 
the cyclic pressure at the pneumatic valving system (blade root), there is a phase lag 
due to length, a pressure attenuation due to duct friction loss, and a pressure rise due 
to the centrifugal pumping effect. The duct pressure is obtained as 



and 

P, » P + P, cos(<|> - <t>) + _ sin(\|> - $) (8) 

dr o ic -ls 

where 4> is the phase lag for the pressure pulse at a radial station defined as 

* - ( r - r J deg (9) 

T a. _ root 

duct 

The term P Q is collective pressure, P lc and P lg are cyclic pressures at the blade root. 
The n,j ct and n UC1 are respectively the duct friction loss coefficient and the centri- 
fugal pumping efficiency, respectively. The V t ^ p is tip speed (ftR), r root is root 
radius (where the pressure duct starts), P m i g atmospheric pressure (1827 lb/ft ) and 
a duct speed of sound in the duct (1274 ft/sec). 

Equations (7) through (9) show that the blowing momentum coefficient is a function 
of radial position r as well as azimuth angle t|». 


VEHICLE TRIM SOLUTION 

Two types of trim solutions are considered corresponding to separate CC rotor 
aircraft concepts. The propulsive trim solution represents a CC rotor implemented in a 
conventional helicopter flight mode. There, the rotor produces all lift and propulsive 
forces. Alternatively, the constrained trim solution represents a CC rotor employed on 
a compound hleicopter configuration that also features auxiliary propulsive devices. 


Propulsive Trim 

The propulsive trim simulates the free flight condition. For specified weight 
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coefficient C^, collective pitch 0 , and forward speed y, the trim solution calculates 

blowing settings (P p and P, ), steady flap response (B , , and 6, ), vehicle 

orientation (a and i ) and steady inflow ratio X* ^ s 

s s 

The present propulsive trim is calculated from the satisfaction of three forces 
(vertical, horizontal and lateral) and two moments (pitch and roll) equilibrium 
equations. Figure 3 shows the forces and moments acting on the vehicle. The 
equilibrium equations have been defined by Panda and Chopra (1985). 

For the vehicle trim solution, only the flap motion up to first harmonic is con- 
sidered. Therefore, the following rotor equations are used. 


(flap equation) di); » 0 (10) 


(flap equation) cosi|> di|> « 0 (11) 
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For steady inflow, a linear distribution model developed by Drees is used. 


C 

X = p tan a_ + — (! + K x £ costJ> + K y f sim|>) (13) 
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For hover, 


K and K 
x y 


become zero. 


As stated earlier, the blowing momentum coefficient varies in both the radial as 
well as azimuth directions. Therefore, the aerodynamic coefficients (C. , C., and C ) 
cannot be expressed in simple analytical expressions. To obtain the rotor coefficients 
(C T , C Y and C^) for trim solution, the element forces must be integrated numerically in 
both the radial and the azimuth directions. Thus, 
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where U , U , U and V are section flow velocity components; see figure 4. The x is the 
nondimens ional coordinate r/R, and o(x) is the local solidity ratio. 

The vehicle and rotor equilibrium equations are obtained for large angles. These 
equations which are expressed in nondimens ional form, are presented in the appendix. 

The nonlinear equations are solved iteratively for the trim solution using the 
Newton-Raphson procedure. 


CONSTRAINED TRIM 

The constrained trim solution, where propulsive force is partially obtained from 
auxiliary power, is calculated by satisfying the vertical force, pitch moment, and roll 
moment equilibrium equations. The vehicle orientation is specified in terms of shaft 
tilt a ; lateral tilt 4> is zero. Since the trim solution is for an isolated rotor, the 
characteristics of the vehicle are not needed. For the specified weight coefficient C w , 
collective pitch 9 , and advance ratio u , the trim solution calculates blowing settings 
(p , P. and P, ) and rotor response (8 , and ®i s )’ This solution procedure is 

similarto that S used in the propulsive trim solution. 

Alternative constrained trim procedure is to fix the collective blowing pressure 
(e.g. P =* 1.5 P ) and adjust the collective geometric pitch to obtain the desired 
thrust level. This trim procedure though not implemented in the present work, would 
yield equivalent blade stability results. 


BLADE RESPONSE SOLUTION 

The blade response solution involves the determination of the time dependent blade 
deflected position. For steady flight conditions, the blade response is periodic, and 
hence the solution is calculated for only one complete cycle. For this, the nonlinear 
coupled blade equations containing periodic terms are solved in the rotating frame using 
an iterative procedure based on Floquet theory (Dugundji and Wendell, 1983). For the 
blade response solution, the pilot controls and the vehicle orientation obtained from 
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the trim solution are used. The numerical procedure to calculate the initial conditions 
and the blade response along the azimuth are detailed by Panda and Chopra (1985). The 
calculated response solution consists of all harmonics for flap, lag, and torsion 
motions . 

It should be noted that there is another cycle of iterations (typically two) between 
the trim solution and the blade response solution to include the effect of blade elastic 
twist in the trim calculations. 


STABILITY SOLUTIONS 

The stability of blade perturbation motion about Its steady deflected position is 
examined using the vehicle trim and the blade response solutions. For perturbation 
motion, unsteady aerodynamic effects are introduced approximately through a dynamic 
inflow modeling. The dynamic inflow is assumed to be a perturbation about the steady 
value, and its components are related to rotor perturbation forces and moments (Panda 
and Chopra, 1985). 

The blade perturbation equations are transformed to the fixed reference frame using 
Fourier coordinate transformation. These equations, which contain selected harmonic 
terms, are solved for stability using Floquet transition matrix theory. 


RESULTS AND DISCUSSION 

Numerical results are calculated for a four-bladed CC rotor with Lock number y * 5, 
solidity ratio a = 0.13, and zero precone. The blade flap, lag and torsion rotating 
frequencies are 2.3/rev, 2.6/rev and 18/rev, respectively. The chordwise offsets of the 
center of mass and the reference aerodynamic center from the elastic axis are considered 
to be zero, and the elastic axis is assumed to be at the midchord position. For stabi- 
lity calculations, the structural damping is assumed to be zero for all modes. For air- 
foil characteristics, tabular data of a typical CC airfoil with trailing edge blowing 
(single slot) is used. Other rotor and vehicle characteristics are given in table I. 


PROPULSIVE TRIM 

Numerical results are calculated for C„/a * 0.1. Figure 5 shows the vehicle pro- 
pulsive trim solution for a collective pitch of zero. The propulsive trim parameters 
P Q , P 1C> P 1S* a s » $s > and * are plotted for different forward speeds (in terms of 
advance ratio y). An advance ratio of 0.6 represents a forward speed of about 300 
ft/sec and a maximum tip speed of 850 ft/sec. The root blowing pressures, collective 
(P ) and cyclic (P^ c and P lg ), are presented in terms of atmospheric pressure (P^). The 
flap angles (8 q> 8 1C and 8^^) are negligible for this highly stiff rotor and hence are 
not presented. The trim solution is calculated iteratively from nonlinear equilibruim 
equations. As conventional rotor, the shaft has to tilt more forward at larger y in 
order to compensate for the increasing parasite drag. The Inflow X first decreases and 
then increases with forward speed due to the combined effect of decreasing induced velo- 
city and increasing disk tilt (a^) at larger y. There is only a slight influence of y 
on side shaft tilt <(>-. The collective and cyclic blowing pressure requirements with 
advance ratio y appear quite similar to the respective geometric pitch requirements of a 
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conventional rotor (Panda and Chopra, 1985). The cyclic blowing pressures (P lc and P ls ) 
are much smaller than the collective blowing pressure (P ). The periodic variation of 
dynamic pressure Is compensated for by the cyclic blowing components. At low forward 
speeds both P. and P. are of equal Importance because of the nonuniform Induced 
velocity; however at higher speeds, P lg becomes larger than P 1( ,. 


TABLE I - CCR Rotor Characteristics In Analysis 


Rotor radius R 


28.5 ft 

Tip speed ftR 


500 ft/sec 

Chord-to-radius ratio, c/R 


0.1 

Airfoil thickness ratio, t/c 


0.15 

Slot height-to-chord ratio, h 


0.002 

Feather inertia-to-f lap inertia 

ratio, I f /l b 

0.0024 

Reference lift curve slope, a 


5.7 

Vertical eg offset from hub, h 


0.2R 

Duct speed of sound a ^ UC £ 


1274 ft/sec 

Duct friction loss coefficient, 

n duct 

0.15 

Centrifugal pumping efficiency, 

n 

pump 

0.57 

Root radius-to-rotor radius ratio, ( r / R ) root 

0.1 


A word of caution: these and subsequent results are calculated using the linear 
inflow model of Drees, which perhaps underestimates the longitudinal Inflow variation at 
low forward speeds (w <0.15). In addition, blowing in the reversed flow regions Is 
assumed to produce no circulatory lift. 


Figure 6 presents the propulsive trim solution for a collective pitch of -10 
degrees. For this pitch setting, a larger collective blowing pressure is needed to 
achieve the prescribed thrust level. With this negative collective pitch, there are 
changes on other trim parameters (A, a and ♦ -); In fact, their values become almost 
double those of zero pitch. Trim results withSut pneumodynamic effects (pressure loss, 
phase lag, and centrifugal pumping) are shown in figure 6(b). There is a small step 
decrease in collective blowing pressure for all forward speeds because of reduced 
losses. Neglecting pneurao dynamics in the calculations also influences cyclic pressures; 

increases . 


decreases and 


■? 
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In figure 7, the time-dependent position of the blade Is presented for one complete 
cycle. The blade Is set at zero collective pitch, and the propulsive trim solution is 
employed. These results are obtained by solving the nonlinear periodic equations using 
an Iterative procedure based on Ploquet theory. For numerical integration, a fourth- 
order Runge-Kutta method and 240 time steps per cycle (At|> * 1.5 degrees) are used. The 
blade is extremely stiff in flap mode (flap frequency - 2.3/rev) and, therefore, a small 
flap response is expected. At a low forward speed (p - 0.2), there is a very little 
flap response. However, at high forward speeds, the flap response becomes greater 
induced by a larger variation in aerodynamic environment along the azimuth. At p *0.4, 
the flap response consists primarily of the second harmonic with a maximum peak-to-peak 
amplitude of about 0.75 degrees. At high u of 0.6, the flap response is again dominated 
by a 2/rev component with a maximum peak-to-peak amplitude of about 2.5 degrees, which 
is perhaps a large flap response for this highly stiff rotor. 

Figure 8 shows the damping of the low frequency cyclic lag mode for different 
advance ratios and collective pitch. Results are obtained from the eigen solution of 
Floquet transition matrix. Again, for the stability solution, 240 time steps per cycle 
are used for time integration. The eigenvalues represent rotor frequencies in the fixed 
reference frame. For this case, the low frequency lag mode is a regressive mode. The 
damping is presented in terms of the real part of the complex eigenvalue, a . Note 
that <*£ ■ <i)£, where is the viscous damping ratio of the lag mode, and^rn^ is the 

frequency of lag mode nondimensionalized with respect to rotational speed. For the case 
of zero collective pitch, the blade is stable but the damping level in lag damping is 
quite low. However, the inclusion of structural damping will increase blade stability. 

In addition, negative collective pitch also has a stabilizing influence on lag mode. 

It is interesting to note that the variation of forward speed has only a slight 
influence on lag damping. 

Figures 9 through 11 show the effect of thrust level on blade lag mode stability. 

For a fixed-collective pitch, the rotor thrust level is a near linear function of 
blowing pressure. Using the previous thrust coefficient of C /a of 0.1. as a reference, 
three additional thrust levels are considered: C /a - 0.05, 0715, and 0.2. Figure 9 
presents the lag mode damping for zero collective pitch. In general, the reduced thrust 
level stabilzes lag mode at low forward speeds and destabiliizes lag mode at high for- 
ward speeds. For C^,/a ■ 0.15, lag mode is unstable for p < 0.15, whereas, for C_/o - 
0.05, lag mode is unstable for p > 0.41. This observed instability is quite weak in 
nature and can be easily stabilized with the inclusion of a small amount of structural 
damping in lag mode. In Figure 10, somewhat similar results are seen for a collective 
pitch of -5 degrees. For a low thrust level condition of C T /o - 0.05, the lag mode be- 
comes unstable for vi > 0.42. At this thrust level, the solution is not obtained for p > 
0.51. The lateral cyclic blowing component becomes larger than the collect iv blowing 
level; therefore, P rQ becomes less than locally on the advancing side of the rotor. 
Again the damping requirements to stabilize the lag mode are not high, and the expected 
levels of internal structural damping would ensure this stability. Figure 11 shows lag 
mode damping for a collective pitch of -10 degrees. With this high negative pitch, the 
lag mode damping generally becomes more stable. The exception is C T /o - 0.2, for which 
the blade lag mode becomes less stable at high forward speeds. 

Figure 12 shows the effect of torsional stiffness on lag mode stability. The 
earlier results are for a rotating torsional frequency of 18/rev. Figure 12(a) presents 
results for zero collective pitch. If the torsion frequency is reduced to 10/rev, there 
is only a slight effect on lag mode damping. A further reduction in torsional stiffness 
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has an appreciable effect on lag mode stability. For a torsional frequency of 5/rev, 
the lag mode becomes unstable and the Instability Increases at high advance ratios. In 
figure 12(b), results are presented for a collective pitch of -10 degrees. Again, 
reducing torsional stiffness decreases lag mode damping. 

The effect of lag stiffness on lag mode damping is presented in figure 13. The 
earlier results are for a rotating lag frequency of 2.6/rev. As compared with conven- 
tional rotors, this is a case of extremely high lag stiffness. For zero collective 
pitch, figure 13(a), if the lag frequency is reduced to the level of a typical stiff- 
inplane hingeless rotor (v - 1.4), the lag damping is only slightly reduced. If the 
lag damping is further reduced to the level of a typical soft-inplane rotor (v » 0.7), 
again the effect on lag damping again is quite negligible. For a negative collective 
pitch of 10 degrees figure 13(b), the effect of reducing lag stiffness to the conven- 
tional rotor value is quite destabilizing. 

Figure 14 presents the effect of flap stiffness on lag mode damping. The earlier 
results are obtained for a rotating flap frequency of 2.3/rev. In comparing with the 
existing rotors, this is an extremely high flap stiffness. For zero collective pitch 
figure 14(a), if the flap frequency is reduced to 1.5/rev (level of ABC Rotor), the lag 
mode becomes slightly more stable, and even more so at high advance ratios. If the flap 
frequency is further reduced to 1.1/rev, to the level of a typical hingeless rotor(e.g. 
BO-105), the lag mode becomes unstable at high advance ratios (y > 0.48). For a collec- 
tive pitch of -10 degrees Figure 14(b), the trends of lag mode damping with decreasing 
flap stiffness are somewhat different. The blade, however, remains stable. 

Figure 15 shows the effect of pneumodynamics on lag mode damping. The blade is set 
at a collective pitch of -10 degrees. Neglecting pneumodynamics effects reduces lag 
mode damping slightly. 


CONSTRAINED TRIM 

Constrained trim calculates rotor controls to achieve a desired thrust and shaft 
orientation. This type of trim condition is possible through an auxiliary propulsive 
device. The solution is obtained by satisfying three rotor equilibrium equations 
(vertical force, pitch moment, and roll moment). This is an isolated rotor trim solu- 
tion, and the airframe characteristics are not needed. With an auxiliary propulsive 
device, it is possible to achieve high forward speeds without causing excessive shaft 
tilts. Thus, the subsequent results using constrained trim include a larger range of 
advance ratios (up to 1.0). An advance ratio of 1.0 represents a forward speed of about 
500 ft/sec and a maximum tip speed of 1000 ft/sec. The compressibility effects, 
however, have not been considered in the present work. 

The constrained trim solution for zero collective pitch and zero shaft angle is 
shown in figure 16. The collective and cyclic components of blowing pressure at blade 
root are quite similar to those obtained with propulsive trim for zero collective pitch; 
Figure 5. The collective pressure P 0 is somewhat smaller at higher advance ratios 
because of reduced inflow through the disk. For y > 0.41, the solution is not 
practical because of the pressure constraint (Proot < Poo) • In figure 17, the con- 
straint trim solution is presented for zero shaft angle and a collective pitch of 
-5 degrees. Once again, the collective pressure requirements at high y are compara- 
tively smaller than those of the propulsive rotor. 
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Figure 18 presents the blade flap response for one complete cycle. For this cycle, 
the rotor is set at zero shaft angle and zero collective pitch. For a low advance ratio 
(y » 0.2), the flap response amplitude is small. At high advance ratios, there is con- 
siderable flap response, consisting primarily of a 2/rev component. For y * 0.5, the 
peak-to-peak amplitude is about 1 degrees, whereas, for y * 1.0, the peak-to-peak ampli- 
tude is about 3.5 degrees. 


Figure 19 shows the lag mode stability results for zero shaft angle. The lag mode 
is less damped for zero collective pitch. The negative collective pitch stabilizes lag 
mode damping. For collective pitch of 0 degrees the results are discontinued for y 
larger than 0.42 because of the blowing pressure constraint ( p root < P^). 

In figure 20, the lag mode damping results are presented for a shaft angle of 5 
degrees. This is a forward tilt of rotor shaft and the inflow through the rotor disk 
Increases, resulting in an increased collective pressure requirement. Compared with 
results obtained for zero shaft angle (Fig. 19), the lag mode is slightly less stable 
for both cases of collective pitch. Note that the zero collective case can now be 
extended up to an advance of 1.0. 

Figure 21 shows lag mode stability results for a shaft angle of -5 degrees. This is 
a rearward tilt of rotor shaft, and the inflow through the rotor disk decreases, 
resulting in a decreased collective pressure requirement. Compared with results 
obtained for zero shaft angle (Fig. 19) the lag mode is more stable for both cases of 
collective pitch. With the reduced collective pressure the range of y for zero collec- 
tive pitch is now reduced to 0.31. 


CONCLUSIONS 

Aeroelastic stability of a simple, three-degree-of-f reedom (flap-lag-torsion) CCR 
blade model in forward flight is examined. Results are obtained using propulsive trim 
as well as auxiliary power trim. Based on the results of this study, the following 
conclusions are drawn. 

Flap response consists primarily of 2/rev, and its amplitude increases with y 
(peak-to-peak amplitude of 2.5 degrees at y * 0.6 for propulsive rotor) 

Lag mode damping becomes more stable as collective pitch is decreased. 

At high forward speeds (y > 0.4), lag mode becomes unstable at low thrust levels 
(C T /cr - .05). 

Reducing the stiffness of a highly stiff CC rotor blade to the level of a typical 
hingeless blade can cause lag mode instability. 

The results obtained with auxiliary power trim are quite similar to those obtained 
with propulsive trim. 
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Appendix 


F(l) 


F(2) 


F(3) 


F(4) 


F(5) 


F( 6 ) 


F(7) 


Nonlinear Equations for Trim Solution 


-(v g -l)/Y 


. , 1 2 f 1 

, B, . - X /h cosct + sina - u 7 — 

h „ /„ 1C eg s s 2 A 


2 £ C T /aa 


w 


< cosa s - -£8. slna g ) + ^ 


(Pitch moment) 


1 - — cos$ cosa - — sina g + — sin<t> g 

w w w 


+ r sln *s + i 1 ‘ 2 ir- ain * s 

w w 


(Vertical force) 


1 2 f 1 

y ■=— + t— cosa - — sina cos* 

2 AC C sC s s 

w w w 


(Longitudinal force) 


Gyp ^y) 

(_ + —> cos<t> g + — sin* g coso g 
w w w 


(Lateral force) 


-(v'-D/y 1 C YP 

2 ic /,; '*• + <c08 * s * •**' > V 

R T 


Y M 

+ sln * 8 - cos* g + ^ (Roll moment) 

C 

X - u tana g - ^ 1 d + f cos* + K y | sin*) (Inflow) 

/y 2 + X 2 

e o * *2- f J 1 x v2 C* dx d * (plap 8 o ) 

v e o o 


327 



2ir 1 2 

F( 8 ) - e ic - - / f x V C t cos* dx d* 


4ir , 2 , . 

(v - 1 ) o o 

8 


2w i 2 

F(9) * B 1<s - 7-7 —5 / / * V C t sin* dx d* 

1S 41T (v - 1 ) o o 

8 


F < 10 > “ 5 1 " f 1 ° ° T V C £ dx d* 

w woo 


C 2x1 

F(ll) - -T - hk- / / « [Up v C t + U T v c d ) sin* 

WOO 


w 


+ (U R V C d - v 2 8 C # ) cos*] dx d* 


C 2ir 1 

F(12) » / * hr - f f 0 [-<°p V C * + U T V V COS * 


C 4if C 
w woo 


+ (U R V C d - v 2 5 C ( ) sin*] dx d* 


(Flap 8 


1C) 


(Flap 8 lg ) 


(Rotor thrust) 


(Rotor drag) 


(Rotor side force) 


In the expressions, the flow components U p , U T and V are nondiraensionalized with 
respect to flR. 
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JET 



QUASI-ELLIPTIC 


Figure 1. - Circulation control airfoil. 


direction^ 
ROTATION 



ROTOR BLADE OPERATES IN FIXED PITCH 

CONTROL AIR EXITS FROM CONTINUOUS 
TRAILING EDGE SLOT 

CYCLIC CONTROL OF AIRFLOW GENERATES 
CYCLIC CONTROL OF LIFT 


Figure 2. - Circulation control rotor concept. 
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COLLECTIVE PITCH =0 DEGREES 

of# = FORWARD SHAFT TILT 

<L = LATERAL SHAFT TILT 
(ADVANCING SIDE UP) 

X = INFLOW RATIO 

(HUB PLANE) «s 



0.2 0.3 0.4 0.5 0.6 

ADVANCE RATIO, M 


(a) Rotor altitude and mean inflow. 

Figure 5. — Vehicle propulsive trim solutions for a collective pitch of zero degree 

(C Jo - 0.1, v„ - 2.3. v_ - 2.6. v- - 18} 


ROOT PRESSURE RATIO 



(b) Blowing pressure requirements for trim. 


Figure 5.- Continued. 
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RADIANS 



(a) Rotor altitude and mean Inflow. 

Figure 6. - Vehicle propulsive trim solutions for a collective pitch of -10 degrees. 

(C T /o - 0.1, Vg * 2.3, v - 2.6, v 0 - 18) 
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ROOT PRESSURE RATIO 



ADVANCE RATIO, M 


(b) Blowing pressure requirements for trim. 
Figure 6.- Concluded. 


334 
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UNSTABLE 

i I I 1 1 — 

0 0.1 0.2 0-3 0.4 0.5 0-6 

ADVANCE RATIO, H 

Figure 8. - Effect of collective pitch on low frequency cyclic lag mode. 

(C T /a » 0.1, Vg - 2.3, v ? * 2.6, v Q * 18) 
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Figure 11. — Effect of thrust level on low frequency cyclic lag mode. 

(v g - 2.3, - v Q - 18) 
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LAG DAMPING 



(b) Collective pitch -10 degrees . 
Figure 12.- Concluded. 
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LAG DAMPING 



ADVANCE RATIO, M 

(b) Collective pitch -10 degrees. 
Figure 13.- Concluded. 


3^3 


LAG DAMPING 



ADVANCE RATIO, M 

(a) Collective pitch = 0 degrees* 

Figure 14. - Effect of flap stiffness on low frequency cyclic lag mode. 

(C^/o - 0.1, = 2.6, Vq — 18) 
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LAG DAMPING 



ADVANCE RATIO, M 


(b) Collective pitch -10 degrees. 
Figure 14.- Concluded. 
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ROOT PRESSURE RATIO 


2.0 



Figure 16. - Blade root blowing pressure requirements for the constrained trim solution. 

(C T /a - 0.1, v g - 2.3, v ? - 2.6, v 0 - 18) 
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ROOT PRESSURE RATIO 



ADVANCE RATIO,/! 


Figure 17. - Blade root blowing pressure requirements for the constrained trim solution. 

(C T /o = 0.1, v g - 2.3, - 2.6, v 0 - 18) 


3H8 


I 



Figure 18. - Blade flap response for various advance ratios. 

= 0.1, Vg = 2.3, = 2.6, v 0 = 18) 


349 



LAG DAMPING. 



0 0.2 0.4 0.6 0.8 1.0 

ADVANCE RATIO, /i 


Figure 19. - Effect of collective pitch on low frequency cyclic lag mode. 

(C T /o - 0.1, v g - 2.3, - 2.6, v Q - 18) 


350 


LAG DAMPING, 



Figure 20. — Effect of collective pitch on low frequency cyclic lag mode. 

(C T /cr - 0.1, - 2.3, - 2.6, v 0 - 18) 
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